Universality and Symmetry Breaking in Conformally Reduced Quantum Gravity 
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The scaling properties of quantum gravity are discussed by employing a class of proper-time 
regulators in the functional flow equation for the conformal factor within the formalism of the 
background field method. Renormalization group trajectories obtained by projecting the flow on 
a flat topology are more stable than those obtained from a projection on a spherical topology. In 
the latter case the ultraviolet flow can be characterized by a Hopf bifurcation with an ultraviolet 
attractive limiting cycle. Although the possibility of determining the infrared flow for an extended 
theory space can be severely hampered due to the conformal factor instability, we present a robust 
numerical approach to study the flow structure around the non-gaussian fixed point as an inverse- 
problem strategy. In particular it is shown the possibility of having a spontaneous breaking of the 
diffeomorphism invariance can be realized with non-local functionals of the volume operator. 
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I. INTRODUCTION 

In statistical mechanics universality is the property for 
which, close to a continuous phase transition, the long 
range behavior of the system is independent of the de- 
tails of the microscopic interactions. For example the 
magnetic order parameter of the Ising model does not 
depend on the lattice geometry, although the critical tem- 
perature is different for a square, triangular or an hexag- 
onal lattice. When different models share the same set of 
critical exponents it is said that they belong to the same 
universality class. 

In quantum gravity, the conceptual difficulty in ex- 
tending the notion of universality is the requirement of 
"background independence" because it implies that the 
geometrical structure of the spacetime cannot play any 
role in the definition of the microscopic degrees of free- 
dom. 

The strategies proposed so far to quantize gravity have 
different ways of dealing with this issue [H-[ll ■ For exam- 
ple, in string theory, although the background indepen- 
dence is not manifest at a perturbative level, it should be 
realized non-perturbatively via AdS/CFT. On the other 
hand, in loop quantum gravity, this requirement is satis- 
fied from the very beginning, at least at a formal level, 
and in the asymptotic safety program [5rl2o| background 
independence is dynamically achieved via the background 
field method [27j . From this point of view conformally 
reduced gravity, a scalar analogous of gravity in which 
the only propagating degree of freedom of the metric is 
the conformal factor, is an important theoretical labora- 
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tory to understand the issue of background independence 
and renormalizability in this context. In fact, although 
in classical general relativity the conformal factor is not 
a propagating degree of freedom, in quantum gravity its 
fluctuations dominate the path-integral because of the 
conformal factor instability and can become the most 
important ones both in the ultraviolet (UV)region, close 
to the non-gaussian fixed point (NGFP) [3, |6(, and in 
the infrared sector (IR), below the gaussian- fixed point 
(GFP). 

The idea of considering this privileged point of view in 
order to better understand the structure of the UV criti- 
cal manifold of quantum gravity was first put forward in 
[HI and . In the latter paper a non-perturbative flow 
equation in the so-called "local potential approximation" 
(LPA) has been derived for the first time and further in- 
vestigated in [30j within the framework of the bimetric 
truncation. Moreover in [3l[ the contribution of the trace 
anomaly and the R 2 term have been considered. 

Important points which deserve to be better investi- 
gated in this approach are the dependence of the critical 
quantities from the threshold functions and the structure 
of the renormalization flow beyond the simple confor- 
mally reduced Einstein-Hilbert (CREH) truncation. In 
the first case it is interesting to investigate the impact of 
the reduced degrees of freedom as a function of the cutoff 
structure. In particular we would like to discuss the pos- 
sibility of determining an "optimal cutoff' for which the 
difference between the calculation of the universal prop- 
erties performed in the full Einstein-Hilbert model and 
in the CREH approximation is minimal. In the second 
case it is important to study the evolution of admissible 
initial data which are not necessarily of the CREH form 
in order to see if non-local contributions to the renormal- 
ized Lagrangian significantly deform the structure of the 
UV critical manifold 1 13. 
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The role of non-polynomial truncations can be signifi- 
cant also in the deep infrared region. For instance in the 
more familiar A ^ 4 -theory in the broken phase, a first- 
order phase transition occurs in three dimensions, where 
the inverse susceptibility is not continuous in the k — > + 
limit. In this case the finite jump of the renormalized 
mass as a function of the field strength is not accessi- 
ble to the standard /3-function approach; moreover only 
with a robust and accurate numerical integration scheme 
of the flow equation it is possible to recover the convexity 
property of the free energy in the thermodynamic limit 

Hill. 

In this work, in particular, we shall use the proper-time 
flow equation which has been extremely successful in the 
calculation of the critical exponents in the 3-dimensional 
Ising model [34| and in quantum gravity [151 ] . The main 
reason to employ this flow equation is that it is rather 
simple to implement different threshold functions and 
to interpolate between a sharp momentum cutoff and a 
sharp-proper time cutoff, as we shall see. 

A question which naturally arises in the context of con- 
formally reduced gravity is if a phase of non-zero mean 
conformal factor (y) 7^ takes place at low energy: this 
symmetry breaking phenomenon could be interpreted 
as a phase of broken diffcomorphism invariance where 
(duv) 7^ and the spacetime geometry naturally emerges 
as a low-energy phase. We shall show that, although the 
IR flow of the LPA cannot be properly determined for a 
continuous set of initial data, it is nevertheless possible to 
study this symmetry breaking phenomenon as an inverse 
problem, within the LPA approximation. In particular 
we shall find a new class of UV fixed potentials which 
evolves towards a low-energy phase where the diffeomor- 
phism invariance is spontaneously broken. 

The structure of the paper is the following: a deriva- 
tion of the proper-time flow functional equation is dis- 
cussed in section II by means of a background "inde- 
pendent" blocking procedure. The fixed points and crit- 
ical exponents are computed in section III for various 
classes of threshold functions and for different projection 
in the background metric. The results are then com- 
pared with those obtained in the full Einstein-Hilbert 
(EH) model. Section IV includes a new numerical dis- 
cretization scheme for the flow equation in the LPA and 
describes the possibility of having a phase of broken dif- 
feomorphism invariance at low energy. Section V is de- 
voted to the conclusions and in the Appendix the explicit 
expressions of the /3-functions for various regularization 
schemes and spacetime dimensions are presented. 



II. WILSONIAN ACTION FOR THE 
CONFORMAL FACTOR 

In this section we shall introduce the concept of Wilso- 
nian action for the conformal factor by means of a con- 
straint average field; this derivation is in fact closer to 
the statistical mechanical point of view than the stan- 



dard approach based on the Schwinger functionals. It 
has the advantage of dealing with a quantity which can 
be directly computed in Monte Carlo simulations since 
the introduction of an external current is not necessary 
(35j in this case. 

Let ^[x] be the action for the fundamental field x( x ) 
which we write as x( x ) — Xb{x) + f(x) where xb(x) is a 
non-dynamical background field and f(x) the dynamical 
(fluctuating) field, and let g^ v be a rigid reference metric 
defined on an Euclidean manifold in d dimensions. The 
Wilsonian action in the presence of the background field 
can be formally defined as 

e -s*[/iXB] = J D [f] S (f k - f) e -^[xB+/] (1) 

where fk(x) is an averaged fluctuation field given as 

/*(*) = J d d y y/jjf(y) p k (y,x; X B) (2) 

being Pk(x,y; y B ) = pk(y,x;xB) a symmetric smearing 
kernel (see [36| for a general discussion on smearing ker- 
nels in Riemannian spaces) whose explicit xb dependence 
needs not to be specified at this level. 

In this formalism x plays the same role of the micro- 
scopic metric 7^ in the full theory. In the complete 
framework a background metric g^ is chosen in order to 
perform the actual calculations, and the fluctuations /i^„ 
are quantized non-perturbatively around this background 
which will be dynamically determined by the requirement 
that the expectation value of the fluctuation field van- 
ishes, {hpu) = hf^, — 0. Any physical length must then 
be proper with respect to the background metric g^ v . In 
the conformally reduced theory the expectation values 
/ = (/) and <f> = (x) — Xb + f are the analogous of 
h^v = (hfu,) and </ M „ = (7^) = g^ u + in the full 
theory. 

The central idea of the conformal field quantization is 
to employ the background metric 

9t»> = Xb %v ( 3 ) 

in constructing the smearing function p^(x,y) = 
Pk{x,y;XB) via the spectrum of — □, being k and k = 
k respectively the momentum operators built with the 
background metric g^ and the fixed metric v = 
2/(d — 2) and d is the space-time dimension. The refer- 
ence metric g^ plays no dynamical role in this process 
but it is fixed to perform the actual calculation, while all 
the dynamical fields are spectrally decomposed using the 
basis of the — □ eigenfunctions whose eigenvalues satisfy 

fc 2 = X B 2 ^ 2 (4) 
in the case of a constant yb ■ The 5- function kernel in (Q]) 

*(/r/) = nW-/») ( 5 ) 

X 
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reduces to S(f — /) in the k —> oo limit, and projects on 
the zero-momentum mode in the k — > + limit instead. 

Therefore at k = the blocked action (fl) coincides 
with the effective potential, namely the non-derivative 
part of the Effective Action, but for k ^ the functional 
described in (p} is an effective action for the "low-energy" 



modes with momentum p < k. The relation with the 
standard Effective Average Action for the conformal fac- 
tor defined in [28|, [29[ can be obtained by noticing that 
the expectation value of the blocking field f(x), which is 
defined as 



D[f(x)} e -s* [/>); xb(*)]; (:e) = J D[ f {x)] D[m] S (f k{x ) _ /») e ^feW+/W] 
D[f(p)]e-*W«)+f{*)\f k ( x ) = (f k ( x )), 

I 



(6) 



can be rewritten, after having introduced the source J{x), 
as 



1 dW k [J; X B] 



Vg dJ{x) 



j=o 



(7) 



where 



W k [J; X B}= d d x{J(x)f(x)-S k [f; X B]}- (8) 



The ijHJ) corresponds to the standard generating func- 
tional W k [J]XB\ defined in [Hj]. 

The (5-function constraint in (fTJ) can be conveniently 
evaluated in the momentum space and the "background- 
blocked" action S k = Sk[f;XB\ can be explicitly com- 
puted in the one-loop approximation. The difference 
ASg = S k+Sk — S k can then be evaluated in the infinitesi- 
mal momentum shell between k and k + 5k, where k is the 
"proper" momentum operator built with the background 
metric g^ v . A functional flow equation is finally obtained 
by taking the Sk — > limit and performing a renormal- 
ization group improvement of the resulting expression 
After this step is accomplished, the "background- 
independent" flow is obtained expressing all the running 
"proper" momenta in terms of the reference energy scale 
k. Rewriting the (regularized) one-loop contribution in 
the Schwinger "proper-time" formalism one finds 



d t S k [f; X B] =~Tr 



— d t T fe exp | - 



S 2 S k [f; XB } 



5p J 

(9) 

where t = log(fc) is the RG time and r k = T k [xB\- 
The important difference between this type of functional 
"proper-time" flow equation and the version used in ear- 
lier investigations [34[ is that the trace in © is here com- 
puted by means of the representation provided by the 
spectrum of — □, 



Tr[A] 



d d xVS (x\A\x) = / d d x^ X B (x\A\x) 



(10) 



The precise relation between the "proper-time" flow 
equation and the exact flow equation 37[ has been ex- 
tensively clarified in 3S] . For actual calculations we shall 
use the one-parameter family of smooth cutoffs r k = r k 
which has been widely used in the literature [34| . whose 
explicit expression reads 



r(n, a Z n k 2 x\ v ) ~ £(n, s Z n Ag utoff x|Q 

r(n) 



(11) 

Here n is an arbitrary real, positive parameter which con- 
trols the shape of the r k in the interpolating regions, 
and r(a, x) = f°°dt t a ~ 1 e~ t denotes the incomplete 
Gamma-function. Furthermore, Z is a constant which 
has to be adjusted to make sure that the eigenvalues of 
— □ are cut off around ~ k 2 rather than ~ k 2 /Z a fl5l |. 
Therefore the derivative dtT k {s) in © explicitly reads 



d t r-(s) lim k 

8k— ¥0 



Sk 



— (Zsk<nxZTexp(-Zsk<nxg) (12) 

with n > d/2. For n = d/2 the kernel (112|) does not regu- 
late completely the UV because the proper-time integral 
requires a field independent (vacuum) contribution to be 
subtracted from the RHS of Eq.([9]). On the other hand 
for n > d/2 this class of regulators allow us to take the 
formal limit n — > oo which coincides with a sharp-cutoff 
introduced in the proper-time cutoff 39]. To conclude 
this section we would like to remark that for n = d/2 + 1 
our regulator corresponds to the "optimized cutoff" in- 
troduced in El. 



III. POLYNOMIAL TRUNCATIONS 

In this section we shall discuss the structure of the 
NGFP obtained by the flow equation © as a function of 
the cutoff parameter n for different reference topologies. 
Let us assume that the field dependence of the blocked 
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action is completely encoded in a relation of the type 
<j) = xb + /j where <j> is a blocked field, so that is a 
local function of <f>. An important example of this ap- 
proximation is the conformally reduced Einstein-Hilbert 
truncation. It is obtained by setting g fLU — 4> 2u in 
the Euclidean Einstein-Hilbert action [4J| 

SfM = -^ Jd d x^gG^{R(g)-2A k ) (13) 
so that the standard formulas for Weyl rescalings yield 



with 



S k [</>] = d d xV¥ Z k \ g^d^ d v <j> + ± A{d) R<f> 2 + 



-2A(d) A fe 0<#W) (14) 



where R = R(g) and 

1 d-1 

Zk. 



2itG k d-T 



A(d) 



d-2 
8(d- 1)' 



(15) 



In particular the second functional derivative of ([14"]) 
reads 



Z k + 2 A(d)R + (16) 

-2 A(d) B(d) A fc <t>T*% " 2 ) y) 



where <J (a;, y) is the Dirac delta function in d dimensions 
in the fixed metric space endowed with the g^u metric, 
and 



B(d) 



2d 



d-2 \d-2 



2d 



- 1 



(17) 



A. 5"* topology 

Let us first consider the topology of the d-dimensional 
sphere S . In this case the curvature of the reference 
metric g^u is constant and the running of the adimen- 
sional coupling g k — G k k d ~ 2 can be obtained from the 
</> 2 term setting Xb{ x ) = constant in the action ([14")) , so 
that 



Si [xb]= Jd d xVf Z k ^A(d)R x i-2A(d)A kX g- 

(18) 

The trace on the background metric can be computed 
using the Seeley-Gilkey-deWitt heat kernel expansion up 
to linear terms in the reference curvature by using the 
Mellin transform Q n , so that 



Tr[W(-D)] = 



(19) 



(4 



^ |Q| J d d x^~g+ l -Q,_ 1 J d d x^~gR^ 



[W{-U)\ 



1 



q n - i W(q)dq. (20) 



r[n] Jo 

In particular, using the relation — □ = — Oxb") 
W(-D) = e - sZk{ -^ = e-^fcX^C-n)), Therefore 



Qi-fJ d ]I dqq " lersZkX ^ q 



i 



{sz k xi u Y- 



(21a) 



|-2 p -»Z hX l v q 



dqq? e 



(21b) 

By inserting expression (jTHJ) and (fT2"j) in the flow equa- 
tion © and using ([T^| with (|2"Tj) . the coefficients of the 

Xg and Xb^ 2 ' t erms due to the renormalized flow of 
the \fgi_iv R (g^ ) and \fg\Z operators are easily identi- 
fied. At last the /3-functions for the adimensional running 
Newton constant g k and the adimensional Cosmological 
constant X k = A k k 2 can be obtained with the introduc- 
tion of the "anomalous dimension" r\ = k d k \nG k , so 
that 

p g {g,\) = kd k9k = [d-2 + n)g k (22) 
and 77 = rj(g k , X k ) in general. In four dimensions we have 



I3 g = g k [2 



9k 



n n T[n - 1] 



(n-2X k ) n - 1 6tti>] 

g k n n £[n - 1] 
(n-2A fe )"- 1 6tti>] 

g k nTjn-2] 



(23a) 



' {n-2X k ) n - 2 2vrr[n] ' ^ 23 ^ 

The expressions for the (i-dimensional /3-functions are 
listed in the Appendix, together with the limiting cases 
n — > d/2 and n — > 00. 

B. R d topology 

In the case of a flat R d topology the scalar curvature 
of the reference metric vanishes, constraining the mass of 
the action (jT5J) to be zero. In order to extract the beta 
function from the flow equation (J5J it is convenient to 
consider a general truncation of the type 



S fc [0] = Jd d x {\g^ Zkd^d^ + V k {(j>)) (24) 

where V k (<j>) = Z k U k ((f>) and employ a derivative- 
expansion around an homogeneous background plus a 
fluctuation, so that <fr — xb + f(x). In this case we have 

S k [<t>] = j d d xVd{-^Z k f(x) □ />) + V k [ XB ] + (25) 
+Vl[ XB ] /» + \ V k n [XB\ f(x) 2 + 0(f(xf) + 0(d 4 /)}. 
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Therefore 

d t S k [f(x)}= (26) 

= -\ jd*x^g X %» J £ a t r» (ale- \ x ) 

where 

K = -Z k U + V k "[ XB ] (27a) 

^ = V fc "' [Xb]/» + \ V k ""lXB} f(xf. (27b) 



The trace in ([26]) can be evaluated in a background- 
independent way by means of an integration in momen- 
tum space over the eigenvalues p 2 of the laplacian built 
from the background metric g^ u , inserting in (|26p the 
identity / d d p\p)(p\ — I(2 7r) d and using in (127a) the 
substitution — □ — > — DXb^- I n or der to disentangle the 
trace in ([26)) a Baker-Campbell-Hausdorff expansion of 
the heat kernel is performed, so that 



dt S k [f{x)]=- 1 - J d d xJl x i» J JUL J ^d t r^s)(x\p)(p\e~ sK (l- S SK+^{[5K,K} + 6K 2 } + ...)\x). 

(28) 



where the dots stand for the higher order terms in the s 
expansion of the exponential and 



-ipx 



(29) 



The matrix elements of the expanded heat kernel can 
then be calculated ordering the operators by means of 
the commutation rule 



(30) 



It is then straightforward to identify the coefficients of the 
V k and Z k terms, obtaining the following set of coupled 
equations 



kd k V k =M(k 2 X B ) 5 



kd k Z h = N{k 2 xl")i~ 



1 



1 4. y fc "(XB) 



w/z k y 



^"(Xb) 

k 2 nZ k x'i 



n+3- 



(31a) 



- (31b) 



where 



M 



n ^# r(n - |) 

r(n) 



(32a) 



JV 



V47T/ 

_ (d - 2 (n + l))(d - 2 (n + 2)) ( n }i T(n - f ) 

(32b) 



24 dn 2 



V47T 



The /3-functions for the adimensional couplings of the 
CREH truncation are then obtained introducing a poly- 
nomial ansatz for the dimensionful potential U k of the 
type 



u k [xi 



k 2 — v 4 



(33) 



so that one obtains in four dimensions the coupled set of 
equations 



f3 g = g k ( 2 - 2 
0x-- 



g k \j I>+l]n" 
(n-2A fc )"- 1 9 7rr[n] 

3fc r(n - 2) 7i" 



+ 



(n-2A fe )"- 2 2 7rr[n] 

3fc A| I>+l]n* 



(34a) 
(34b) 



+A fc -2-2 



(n-2A fe )"- 1 9 7rr[n] 



Fixed Points and linearized flow 



The /3-functions (J23J) and flMJ) vanish both at the GFP 
located at A* = = 0, and at a NGFP defined at 
A* 7^ 0, g* 7^ 0. The properties of the linearized flow 
around the NGFP are determined by the stability matrix 
B 



Bij=[d gi /3 



y, 



{ff>} = {9*} 



(35) 



{9i} G {g, whose eigenvalues 6*i i2 = —8'±i9" form in 
general a complex conjugate pair. A negative real part of 
the eigenvalues, i.e. a positive 9' (we will refer to it as the 
first Ljapunov exponent, following the standard notation 
used in dynamical systems), implies the stability of the 
fixed point, while the imaginary part characterizes the 
spiral shape near the fixed point. Our results in four 
dimensions are summarized in Table HI while results in d 
dimensions can be found in Table [TT] in the Appendix. 

It is clear from Table U that also the theory defined 
by the CREH approximation is asymptotically safe, al- 
though the scaling properties are rather different from 
those obtained from the full EH in [13 . For instance the 
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TABLE I: The fixed point values and the critical exponents obtained in four dimensions from the /3-functions (|23a|) an d (I23b[) 
for various values of the cutoff parameter n compared with those previously obtained from the full EH gravity in [13], on the 
right. 



CREH-S 4 



CREH- 



Full EH-S" 4 



A* <y* 



0' 



\*g# 



0' 



A*<7# 



0" 



3 
4 
5 

6 

8 
10 
15 
20 
30 
50 
100 
300 
00 



1.125 
1.2 
1.25 
1.285 
1.333 
1.364 
1.406 
1.428 
1.451 
1.470 
1.485 
1.495 
1.5 



1.571 
1.810 
1.885 
1.914 
1.930 
1.931 
1.923 
1.914 
1.904 
1.894 
1.886 
1.880 
1.880 



1.767 
2.171 
2.356 
2.461 
2.574 
2.633 
2.703 
2.735 
2.764 
2.785 
2.800 
2.810 
2.815 



3 

1.5 
1 

0.75 
0.5 
0.375 
0.230 
0.167 
0.107 
0.062 
0.030 
0.010 




4.795 
4.213 
3.873 
3.665 
3.427 
3.295 
3.129 
3.050 
2.974 
2.914 
2.871 
2.842 
2.820 



0.800 
0.837 
0.867 
0.889 
0.921 
0.941 
0.971 
0.987 
1.004 
1.018 
1.028 
1.036 
1.040 



2.084 
2.666 
2.914 
3.041 
3.159 
3.209 
3.253 
3.265 
3.271 
3.271 
3.269 
3.266 
3.265 



1.670 
2.234 
2.528 
2.706 
2.910 
3.023 
3.161 
3.225 
3.285 
3.331 
3.364 
3.385 
3.396 



8.580 
5.721 
5.000 
4.578 
4.102 
3.839 
3.511 
3.356 
3.206 
3.089 
3.003 
2.948 
2.920 





2.928 
3.428 
3.627 
3.788 
3.850 
3.903 
3.919 
3.929 
3.933 
3.934 
3.934 
3.923 



0.355 
0.265 
0.230 
0.211 
0.191 
0.181 
0.169 
0.163 
0.158 
0.154 
0.152 
0.150 
0.103 



0.388 
0.472 
0.517 
0.546 
0.582 
0.603 
0.630 
0.644 
0.658 
0.668 
0.676 
0.682 
0.685 



0.138 
0.125 
0.119 
0.115 
0.111 
0.109 
0.106 
0.105 
0.104 
0.103 
0.102 
0.102 
0.070 



1.835 
1.770 
1.750 
1.742 
1.734 
1.731 
1.727 
1.725 
1.723 
1.722 
1.721 
1.721 
1.720 



1.300 
1.081 
1.000 
0.959 
0.916 
0.894 
0.868 
0.856 
0.846 
0.837 
0.831 
0.828 
0.826 



critical exponents 9' and 9" display an n-dependcncc 
which is stronger in the case of the CREH than for the 
non-reduced theory, although the quantity A* is rather 
stable in both cases. 



We can quantify the impact of the EH conformal re- 
duction with respect to the full EH theory by defining 
a x 2_ type of "distance" in the space of the "universal" 
quantities, by means of 



J 



X 2 (n) = 



(Kg*(C) - Kg*(E)) 2 , (9'(C) - 9'{E)f , (9"(C) - 9"{E)f 



9'{cy 
I — 



-9'{E) 2 9"{C) 2 +9"{Ef 



(36) 



where "C" and "E" stands for CREH and EH, respec- 
tively. 

A plot of this quantity as a function of n is depicted 
in the upper panel of Fig.([T]), for the S 4 projection (solid 
line) and the K 4 projection (dashed line) where it is clear 
that the minimum is attained for n = 4 in both cases, 
which corresponds 1 to the optimized cutoff introduced in 
[40| . On the other hand, in the case of the W 1 projection 
the scaling properties are much less sensitive to the cutoff 
parameter n and the n = 00 limit is as good as the n = 4 
case. 

Of particular interest is the n = 00 limit for the S d 
topology, in which the first Ljapunov exponent vanishes. 
In this case the theory is still UV finite although not 
asymptotically safe anymore, since now the linearized 
system is defined by pure imaginary eigenvalues ±9" and 



1 For an extensive review about the correspondence between ERG 
and proper time cutoff functions see [3H . 



every perturbation of the NGFP will evolve in a cyclic 
trajectory. 

It is also interesting to discuss the scaling properties 
of the theory in the S d projection as the dimension is 
changed. This is shown in the central panel of Fig. (JT|) 
for n = 4 for 9', 9" and for the dimensionless quantity 



Td 



2/(d-2) 



A fe G 



2/(d-2) 



The first Ljapunov 



exponent 9' vanishes for a critical dimension value d c so 
that the fixed point undergoes an Hopf bifurcation as the 
dimension d crosses d c (represented in Fig. ([2])). 

As it is shown in Fig.©, for d — >• 2 the cycle col- 
lapses on the A = line. In this regime it shows a non 
homogeneous running due to the low transient of the tra- 
jectory near the GFP, while it becomes an homogenous 
slow transient around the NGFP in the limit d — > d c . 

Notice that the critical dimension is a function of n, 
d c = d c (n), and while for n — 00 the critical dimension 
is d c = 4, generally holds d c (n) < 4 for a finite value of 
the parameter n. At d = d c the UV behavior is regulated 
by a limit cycle whose behavior resembles the one of the 
Van der Pol oscillator [441 ]. 
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2 3 




FIG. 1: Upper panel: the quantity X 2 ( n ) as a function of 
the cutoff parameter n in the case of S 4 projection (solid line) 
and R 4 dashed line. Central and bottom panel: the quantity 
Td (solid line) , 0' (dashed) and 9" (dot-dashed) as a function 
of the dimension d for n — 4 in the case of 5 4 (central) and 
M 4 (bottom) projections. 



the trajectories which survive for d > d c and that require 
higher-order operators in order to cure the IR sector. The 
second region is the set of points inside the cycle which 
flow towards it in the UV and towards the NGFP in the 
IR. The latter case leads to a new interesting scenario in 
which the UV and IR critical manifolds coincide and the 
EH truncation is finite at every energy scale. 

For this scenario to be plausible we require the cyclic 
trajectory to be close enough to the GFP, so that it shows 
a semiclassical regime. Unfortunately, as it can be seen 
from Fig. ([3]), in the best case (d c = 4 for n = oo) a limit 
cycle with a good semiclassical regime occurs only for 
d ps 3. It is also important to stress that the limit cycle 
never approaches the singularity A = n/2, where the EH 
truncation stops to work. 

Since the Hopf bifurcation is not present in the W 1 pro- 
jection for the CREH, also for small values of the dimen- 
sion, we analyzed the behavior of the linearized flow near 
the NGFP in the case of the full EH truncation, to verify 
if the Hopf bifurcation is still present in the S d projection 
for some value of the parameter n. The /3-functions are 
listed in the Appendix together with results collected in 
Table UH and are a simple d-dimensional generalization 
of the results reported in [l5| . As it can be seen from 
Table |H] the full theory presents a stable NGFP in the 
whole n — d plan, which means that the contribution of 
spin-2 degrees of freedom lower the value of the critical 
dimension under the "critical" value d = 2. 

Although such a non-trivial behavior in the UV region 
seems to be a direct consequence of the strong depen- 
dence of the flow in the S d projection from the cutoff 
parameter n, it is interesting to notice that recent inves- 
tigations based on "tetrad only" theory spaces also show 
the presence of limit cycles in the UV limit 45]. We 
would also like to stress that the nature of the UV limit 
cycle and its Hopf bifurcation discussed in this work is 
indeed different from that of the IR limit cycle recently 
found in the minisuperspace approximation for the EH 
truncation [46] , since in this latter case the IR limit cycle 
originates from the competing fixed points structure of 
the phase-space [47| . 

Finally, it is interesting to note that in the full-EH 
truncation the NGFP is defined for a vanishing cosmo- 
logical constant in d = 3 for every value of the cutoff 
parameter n. 



IV. NON-POLYNOMIAL TRUNCATIONS AND 
SYMMETRY BREAKING 



For d < d c (see left panel of Fig.©) the theory space is 
now divided in two regions. The first is the set of points in 
parameter space outside the cycle, which trajectories flow 
towards the UV to the limit cycle and hit in the IR the 
singularity A = n/2 (or flow towards A = oo). Those are 



In this section we would like to understand the struc- 
ture of the UV critical manifold beyond the polynomial 
truncation discussed in the previous section. We are in- 
terested in the possibility of having a transition to a phase 
of broken diffeomorphism invariance at low energy (29j 
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FIG. 2: Flows in d-dimension for the CREH S projection. Left panel: the limit cycle at d < d c , the dashed line is the repulsive 
internal flow. Center panel: the limit cycle at the critical dimension d — d c . Right panel: the flow of the UV attractive NGFP 
at d > d c . The red line is the location of the FP as a function of the dimension d. The plots have been obtained setting n — oo 
and, starting from the left, for d = 3.9, d — 4, d — 4.1. 



3.5 
3.0 
2.5 
2.0 
1.5 
1.0 
0.5 



the potential by denning the quantity 



FIG. 3: The Hopf bifurcation as a function of the dimension 
d for n = oo in the g-X plane. The solid line is a cycle for an 
initial value near the FP at d c — 4, the other cycles plotted are 
at d — 3.9 (dotted), d — 3.6 (dot dashed), d — 3.3 (dashed) 
and d = 2.9 (long dashed). The red line is the location of the 
FP as a function of the dimension d. 



and thus we intend to numerically solve pip . In fact 
although in four dimensions, an ansatz of the type 

V[ X ] =c 1 (k) X 2 +c 2 {k) X \ Z = c 3 (k) (37) 

is an exact polynomial truncation of the coupled equa- 
tions (|31"j) ; the question concerning the RG evolution of 
a generic initial data determined by (|31D cannot be an- 
swered with this strategy. The core of the problem can 
be grasped in the LPA approximation which solves only 
pial) by assuming a RG evolution for the wave-function 
renormalization functional Z^. In fact although the solu- 
tion of the coupled problem (|3"Tj) is beyond the aim of this 
work, we present a successful numerical strategy to deal 
with (I31a|) which we hope can eventually be extended to 
treat the coupled system pial) and plbl) beyond the sim- 
ple CREH truncation. In particular, in this section we 
shall investigate the role played by higher powers of vol- 
ume operators of the type ~f = (J d d x^Jg) in providing 
a transition to a phase of broken diffeomorphism invari- 
ance. In order to carry out the numerical integration of 
(I31ap it is useful to "linearize" the evolution equation for 



W[ X ] = x 4 1 + 



V" 



nk 2 Z x 2 



(38) 



with 7 = n— 2 > which diverges at +oo as the "spinodal 
line" nk 2 ZkXo v + ^fc'(xo) = is approached, but it 
behaves as a power law for large values of the field outside 
the "coexistence" region where nk 2 Z Xo + V"[xo] < 0. 
In terms of this new variable Eq. (I31ap reads 



(39) 



(2 + rtink 2 Z X 2 (W~~ x~ ~ 1) 
— n k 2 Z x 



k d k W = A n fc 4 di x W. 



The advantage of this manipulation is that equation p9p 
is now linear in the second derivative. 

Ideally we would like to evolve an initial data defined 
at the cutoff scale along the RG direction, i.e. towards 
the infrared. This is usually achieved by defining the 
RG time via k = e _t with t > so that the Cauchy 
problem is fully determined when W[x, 0] is fixed and 
W^IXout, t] is given, Xout being an asymptotic value of the 
field (xout 3> 1 in actual calculations). However if we 
intend to do so, we immediately run into the difficulty 
that as Z < 0, equation (|39|) belongs to the restricted 
elite of the backward-parabolic equations, i.e. a class of 
diffusion-type partial differential equation with a negative 
diffusion constant. As it is well known, in this case the 
Cauchy problem is not "well-posed" (in the sense of [Hj]) 
and the existence of the solution for generic initial data 
is not guaranteed even for an infinitesimal time step! 

Although we already know that an admissible initial 
data in four dimensions is precisely of the type ([3Tj) it 
seems that the CREH truncation is one of the few initial 
data compatible with the general flow equation p~Tj) . In 
fact, the important question concerning the possibility of 
developing a non-zero vacuum expectation value of the 
conformal factor in the k — ¥ + limit cannot be answered 
within the simple CREH truncation. 

It is therefore necessary to treat the question as a sort 
of inverse problem and to consider an integration in the 
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UV direction instead, so that the RG time is defined as 
k = e* with t > 0, the Cauchy problem is well posed and 
the solution is unique. Clearly, once the solution in the 
deep UV is found it is possible to argue that precisely that 
solution is an admissible initial data for a non-singular IR 
flow. However, also in the case of the UV evolution, due 



to its strong nonlinearities, a proper numerical strategy is 
to implement a fully implicit predictor-corrector scheme 
on an uniform spatial and temporal grid. 

In particular, the predictor step is computed at times 
t = (j + 1/2) At so that we can discrctizc (l39l) according 
to the scheme 



^S 2 x W iJ+1/2 = (40) 
^(2 + r,)e- 2 W ( W -/m$ - l) - ^e^( t h)^W^- 12 - (W iJ+1/2 -W id ) 

being h — 1/Ax the spatial grid spacing, q — 1/Ai the temporal grid spacing and, as usual, S 2 Wi = Wi-i — 2Wi + 
Wi+i. The corrector step is instead given by 



2h 2 



€ [W itj+1 + Wi. 



(41) 



n Z, 



»,i+l/2 



(2 + T 1 )e~ 2t {ihf 



IF 



/2 (ih)i - 



1 - 



n Z; 



A 



+2 



and the solution at j + 1/2 in (|4Tj) is obtained from (j40|) 
from the solution of the linear tridiagonal system problem 
in the predictor step. As a consequence (|4"Tj) also reduces 
to a linear problem for the j + 1 time step which can 
be conveniently solved by standard tridiagonal solvers. 
The method is thus unconditionally stable and 0[h 2 + q 2 ] 
accurate (43[. The boundary condition at x = Xinit is of 
the von Neumann type so that d x W = at the inner 
boundary, and < Xinit -C 1: we have checked that our 
results are rather insensitive to the choice of Xinit which 
could be set arbitrarily close to zero in all calculations 
(note however that, strictly speaking, x > always). 

The outer boundary is taken at some x ou t 3> 1 where 
for W a power-law behavior is assumed, like in the more 
familiar Ising model [32| . The initial value of the poten- 
tial at the cutoff reads 

V{xA = ^X 4 + vX 6 +u X 8 (42) 
6 

where the bare values of A, a and ui have been chosen in 
order to display a non-zero minimum as an initial con- 
dition. In addition we have considered also the coupling 
lu < in order to have a real function W for large val- 
ues of the field. In fact, unless lu = a = no consistent 
initial condition can be given in all the real line for the 
potential, as the threshold functions (the denominator in 
pial) ) become complex at a finite value of x for to > 0. 
In solving (|39p close to the NGFP, we have set 77 w —2 
and Z ss — k 2 /g* in d = 4 as we are interested in the UV 
evolution. 



Our results are then summarized in Fig. ((4]): in the left 
panel a symmetry breaking initial state evolves towards 
a convex potential as the UV evolution is followed. The 
final, fixed point state, is then reached already for t = 6 as 
it can be seen in the right panel. Note the "flat bottom" 
of the potential and the almost exponential suppression 
at large values of the field in the final solution. We found 
that the appearance of a "fixed point" potential of the 
type showed in Fig.Q seems to be quite generic if the 
initial condition is changed. 

In Fig. ([5]) another example of the evolution is shown 
for n = 5 and for a different set of initial conditions. Note 
however that the large field behavior of the potential is 
characterized by an inverse power behavior for large value 
of the field thus signaling the presence of non-local invari- 
ants in the fixed point potential. The conclusion of our 
numerical experiment seems to suggest that there exists 
a more complex fixed point structure, not necessarily of 
the of the CREH type, whose precise structure is unac- 
cessible with the more standard /3-function approach. 



V. CONCLUSIONS 

In this work we further explored the universal proper- 
ties of the CREH theory around the NGFP. 

The NGFP characterizes the UV evolution for differ- 
ent projection in the theory space and for a large class 
of threshold functions although we also found the possi- 
bility that the continuum limit is defined by means of a 
limiting-cycle in some cases. Moreover it is possible to 
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FIG. 4: The dimensionless initial potential as a function of the RG time t for an initial condition with A = —1/2, a — 0.05 and 
uj = —0.00714. In the left panel the early RG evolution is showed: t = (dashed line), t = 0.02, (dot-dashed line) and t — 0.8 
(solid line). In the right panel the deep UV fixed point is presented: t — 6 (solid line) and t = 8 (dashed line). Further increase 
in t did not show significant changes in V . 



-0.004 





FIG. 5: The dimensionless initial potential as a function of the RG time t is shown in the left panel for an initial condition 
with A = —6, a = uj — — 1 and n = 5: t = 8 (solid line), t = 32 (dashed line). In the right panel the corresponding evolution of 
the function W is shown. 



find an "optimal" threshold function for which the renor- 
malization flow minimizes the differences with the calcu- 
lation of the universal quantities in full Einstein-Hilbert 
truncation. 

However, going beyond the CREH truncation has 
proven to be rather difficult because of the flawed struc- 
ture of the flow equation in the presence of the confor- 
mal factor instability, and only an UV integration was 
possible. How can we make contact with "real" gravity 
which is formally defined at k — in more general trun- 
cations? Clearly, a mechanism to stabilize the conformal 
factor modes is needed and, most probably, the inclusion 
of the R 2 term is essential to provide a consistent trun- 
cation [48| [49|. Recent investigations discussed in [Ho[ 
for a general f(R) can provide important information for 
a more complete understanding of the renormalization 
flow in the theory space. The structure of the UV region 
around the NGFP can be richer than expected, with a 
class of fixed point potentials displaying an inverse power 



behavior for large value of xb , suggesting the presence of 
non-local volume invariants in the fixed point potential 
14]. Within this class of local potentials it seems pos- 
sible to realize a situation where the diffeomorphism in- 
variance is broken at low energy. It would be interesting 
to further extend this investigation by including a con- 
sistent running of Z and the contribution from the R 2 
term to see the infrared flow structure of a i? 2 -stabilized 
conformal factor beyond polynomial truncations, and we 
hope to address this point in a following work. 
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TABLE II: The fixed point values and the critical exponents obtained in for the cutoff parameter n — 4 from the /3-functions 
(I23al and (|23h>|) for various values of the dimension compared with those obtained from the full EH gravity, on the right. 
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Appendix A: /3-functions 



In this Appendix are listed the explicit expressions of the /3-functions in d dimensions for different choices of the 
cutoff function. The /3-functions obtained using the cutoff (|11[) will be referred as Smooth Proper Time Cutoff, where 
the limits n — ¥ d/2 and n — > oo will be called respectively Sharp Momentum Cutoff and Sharp Proper Time Cutoff. 



1. The projection on S d 



• CREH - Smooth Proper Time Cutoff 



->2-d 
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d-2- 



2 (d-1) 



9k n" 
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2 3 ~ d It 1 -': 
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2 2 - d (d-2)7r 1 -2 5fc n™ri 
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-4 + 1 
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• CREH - Sharp Proper Time Cutoff 
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9k 
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CREH - Sharp Momentum Cutoff 



Pg = 9k 



d-2 + 



d Ti 



(A5) 



/3a = 



2 3 -¥ (d-2)^d d / 2 



.9/,: 



\ 2 3 -¥ d d / 2 Tr 1 -^ <? fc log [l - ( ^"i )A 
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• FULL EH - Smooth Proper Time Cutoff 
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FULL EH - Sharp Proper Time Cutoff 
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• FULL EH - Sharp Momentum Cutoff 
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2. The projection on 
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CREH - Sharp Proper Time Cutoff 
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1 2 3(d-2)3(d-l)3 cxp 



2 (d - 1) 



Afc (A16) 



• CREH - Sharp Momentum Cutoff 



2 5 -¥ (d + 2)3 ^-i d d ' 2 g k X 3 k (3 (d 2 - 3 d + 2) 2 + (d + 2) 2 A 2 - 3 (d - 2) (d - 1) (d + 2) A fe ) 
:// ! " ~ 2 3(d - 2) 3 (d - l) 3 (d 2 - d(A + 3) - 2A + 2) 3 r (f + l) 

(A17) 
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2 5 -¥ (d + 2) 3 tt 1 ^ d d ' 2 g k \ 3 k (z (d 2 - 3 d + 2) 2 + (d + 2) 2 \\ — Z(d—2){d— 1) (d + 2) A fc ) ' 
3 (d - 2)3 (d - 1)3 (d 2 - d (A fe + 3) - 2A fe + 2) 3 r (f + l) 



+ 



2 3 -¥^-! rf ^ ft log 



(d+2)\ k 



(d-2) (d-1) 



(A18) 
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